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ABSTRACT 



A mathematical model is developed for determining an optimum 
submarine tender load list for the POLARIS (FBM) weapon system. 
Military essentiality of spare parts and the stowage capacity of the 
tender play prominent roles in the development of the model. The 
measure of effectiveness which is maximized is the average number 
of demands filled over a given period of time. This is equivalent to 
minimizing the number of stockouts of an item. An example of the 
use of the model is included for the case where the demand for 
items is characterized by the normal distribution. The paper in- 
cludes a computer flow chart and program to aid in implementing 
the model. 

The authors wish to express their appreciation for the 
assistance and encouragement given them by Associate Professor 
Julius H. Gandelman of the U. S. Naval Postgraduate School. 

The authors are also indebted to Miss Patricia Hoang for her assis- 
tance in programming the model on the CDC 1604 computer. 
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z i 

x i 

*i( x i) 

eM 

b(zj) 

s(zj) 

X 

c i 

c 

c 

Mi 

mi 

P 



a' given item where i =1,2,3 N (the entire 

list of spare parts is drawn up and then a number 
is assigned to each spare part to identify it). 

amount of item i to be stocked. 

tVi 

demand for the i l item. 

probability that an amount Xj will be demanded 
during the period under consideration. 

expected demand for item i ; = mi. 

expected shortage of item i. 

th 

average number of 1 n item supplied to meet demand. 

Lagragian multiplier. 

cube (volume) of the i ^ item. 

total cube available for spare parts stowage on tender, 
total cube associated with any A for which 

th 

essentiality (military worth) number assigned to i 
item. 

th 

deviation of i item from the mean, 
mean of the i^ item. 

th 

characteristic number associated with the i item. 
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1. Introduction. 



The trend toward the development and utilization of highly com- 
plex, sophisticated weapon systems within the Navy has generated 
an ever-increasing requirement for more efficient logistical support. 
This paper is an attempt to develop a mathematical model that can 
be used as a decision aid for one aspect of the logistics problem 
associated with the POLARIS system. The model provides a method 
for determining the items and associated quantities that should be 
stocked aboard the supporting tender. 

The requirement for more effective tender load lists is even 
more evident when viewed in conjunction with the POLARIS mobile 
base concept. 

In the event of enemy attack, the Navy's immediate POLARIS 
retaliatory capability will have been dependent on the ability of the 
supporting tender to provide the necessary SSB(N) replenishment 
and repair requirements. This aspect of the tender function is fur- 
ther complicated by the limitations of stowage space available for 
stocking the required spare parts. As a result of this space 
limitation, tender load lists must reflect SSB(N) requirements, and 
tender re- supply procedures must be closely coordinated and 
evaluated to insure that SSB(N) requirements are being met in an 
effective and efficient manner. 
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The concept of military worth is incorporated in the model to 
reflect the assignment of priorities between the various items. 

These essentiality numbers are an attempt to define the military 
worth of a particular item with reference to the performance of the 
integrated weapon system. Essentiality factors will be assigned 
arbitrarily; however, it should be apparent that actual MEC numbers 
could have been utilized for illustrative purposes. 

The measure of effectiveness to be maximized will be the average 
number of demands for spare parts that are filled by the tender over 
some period, as a function of the military essentiality of the various 
items. This is equivalent to minimizing the number of stockouts 
of an item. 

It is obvious that if the demand for an item was known exactly, 
the problem would be simplified considerably subject to the space 
limitation of the tender. Although demand cannot be predicted with 
certainty it may be possible to estimate the probability distribution 
of demand for each item on the basis of past usage data. In the event 
that usage data are not available, the model can be utilized by as- 
suming a probability distribution of demand for the item(s), based 
upon experience and judgment. This would be the situation encoun- 
tered in developing an initial tender load list for a new class of 
submarines. As usage data become available the tender load list 
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can be continually updated to reflect actual demand, utilizing 
automatic data processing equipment. Appendix III contains a com- 
puter flow diagram and program that will accomplish this objective. 

Although this paper is primarily concerned with presenting a 
method for determining an optimal tender load list, certain aspects 
of the tender re- supply problem are included for consideration. 
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2 . 



Assumptions 



The following simplifying assumptions are introduced in order 
to illustrate the proposed method. 

1) The model is concerned with spare parts and does not 
consider dry stores and other consumable items. As a result, C 
refers to the cube available for spare parts stowage. 

2) is a continuous variable, and its associated probability 
density function f^(x^) possesses the properties required for the 
existence of derivatives and integrals. 

3) The demand for each item is independent of the demand 
for every other item. 

4) The probability distribution of demand is known for each 
item over the time period considered. 

5) The list of items considered reflects only those for which 

demand has been generated. If demand for an item stops , then the 
item is removed from stock and a new load list is generated. Con- 
versely, if demand for a new item is generated, it is added to the 
load list and a revised Z = |z^ z 2)' • Z N | determined by means 

of the model. 

6) The development of the model is independent of the re- 
supply phase of the tender operation. The load list determined is a 
function of demand over some period during which the tender must 
provide support without being re- supplied. Theoretically the tender 
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is completely depleted of spare parts at the end of the period con- 
sidered. The effect of re- supply as a function of the desired supply 
objective will be considered in a later section. 

7) The utility values defined in Appendix I are additive. 
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3. Presentation of the Model. 

This section presents the essential concepts of the proposed 
model. Appendix I contains a detailed mathematical formulation of 
the results of this section. 

The expectation, ECU, of x. with respect to the 

probability function PL.*3 gives the following expression for the 
average amount of item i demanded: 

(!) EC*il = f i (*i) dx i» 

where X i= demand for the i^ 1 item. 

(The limits of integration range from 0 to OO since only items 
for which a demand is generated are considered) 

Consider the following situations: 



( a ) X i < z i , 

(b) x i > z^, 



where zj = amount of item i to be stocked. 

If demand for an item is less than the amount of the item that 
isstocked then the tender can meet demand. However, if demand 
exceeds the stock level then only an amount can be supplied. If 
this condition occurs (x^ - z^) represents the shortage (backorder) 

associated with that item. The expected shortage of an item can be 
expressed as a function of the amount of the item stocked. 

(2) b^z^ = J (xi— z^ fj(x i ) dx^ For x^zp 



0 



otherwise . 



The expected value of the amount actually supplied to meet 

demand is given as a function of the amount stocked on the tender. 
This results in the following expression: 

(3) Si £ Zi 3 = J* J 34 fi(x.) dx { + Zi fiUiJdxi 
Manipulation of equations (1-3) confirms the intuitive concept 

that the expected demand must equal the expected demands filled 

plus the expected shortages. * 



(4) ECxil = + bi ^3 

As mentioned previously, the measure of effectiveness is a 
function of the military essentiality of the item as well as the de- 
mands filled by the tender. This requirement is met by defining 

N 

M(z, , . . z n )= 2M. s^(z-) in accordance with assumption (7). 

The problem is now one of determining the set ^=J| •*' 2 .^ 

such that 









is maximized subject to the constraint cf c- z- = C , z^>0. 



Equation (5)(below) is derived in Appendix I and is a general 
result of the considerations mentioned in this section. This equation 
provides the means for determining an optimum tender load list. 

(5) f fi(xi) dxi =l-\(2L), i= 1, 2 N, 

*0 

where A is a lagrangian multiplier with 

M * 

values in the range min ( — ) > A — O 
subject to: ^ 

Jr = C n ^i>o. 

^Appendix I 
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4. General procedure for utilizing the model. 

The general procedure for solving equation (5) is to assume a 

X *A ' 

corresponding to min ( ) , and compute the set Z 

by using the appropriate cumulative distribution function table. The 
C 1 obtained (where C' results from using a A other than that for 
which when the are substituted into the constraint 

equation ( 21 - C- ) determines the minimum required stowage 

space aboard the tender as a function of demand over the demand period 
associated with f^(x-). (As X increases, C 1 decreases. Since 
min ( ) is the largest permissible value of X > the corresponding 

C/is a minimum,) It should be noted that for any \ in the range 
min X — O , every ^ 0. This automatically sets a mini- 

mum stowage requirement on the tender of at least the sum of the 
cube of each item. 

Under the assumption that C is known for the tender, a com- 
parison is made between this value and C*. If C > C 1 an iterative 
solution is required to determine the X which results in C 1 = C. 

The corresponding z^’s obtained are the optimum load list. 

If, on the other hand, C ^ C 1 the model indicates that the ten- 
der cannot effectively meet demands without re-supply during the 
period under consideration. In the event this occurs several decision 
criteria can be adopted. If demand for each item is relatively uni- 
form, one solution would consist of decreasing the tender load list 
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time requirement. For example: Let C 1 = 15,000 cu. ft. represent 
a ninety-day load list minimum requirement. If C = 5,000 cu. ft. , 
then the tender can in effect maintain a 30 day load list by decreasing 
each proportionally. 

This criterion has the advantage of maintaining a load list which 
consists of all items that are demanded. The most obvious disadvan- 
tage is the increased re- supply requirement. 

Another possibility is the elimination of items from the load 
list. This requires a decision criterion to select the item(s) which 

should be deleted. Since C 1 decreases as \ increases .the model 

M • 

indicates that the item corresponding to min( ) should be the first 

item removed from stock. If (C 1 — c^ a^ = C ; a^ < z^) then an 

amount (z^ - a^) can be maintained for this item. If(C ! — cjzi>C) 

then the item corresponding to min ( — ~ ) is eliminated completely. 

M * 

A new min( - gr- ) is selected from the remaining list of N-l items 

and the process is repeated. When a X can be found such that 0=0* 

the z^'s have been determined for this decision criterion. This 

particular policy insures that the items with the highest essentiality 

to cube ratio (__±_) are retained on board. 

C-L 

hhat period during which the tender must be able to operate 
without benefit of re- supply and still maintain effective support 
capability. 

2 Note that other policies could have been chosen. For example, 
one might choose a policy that insures that items with high values of 
P are retained on board, while items with low values of P are eliminated. 
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It permits tender support to be maintained at some reduced effective- 
ness level. (This reduced level can be a significant one). In effect, 
this is analogous to an unexpected demand situation and will be con- 
sidered further in connection with tender re-supply requirements. 

The iterative nature of the solution and the number of items 
comprising load lists indicates the necessity for utilizing electronic 
computers if the system is to be incorporated effectively. Appendix 
III contains a computer description and flow chart for the model. 
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5. 



Example and General Results: 



A numerical example will serve to illustrate the procedure 
for determining an optimal tender load list utilizing the proposed 
model. Consider the case where the demand for each item is 
characterized by the normal distribution. From equation (5) 



-w 



1 



" e.” ' 4 *‘dL A - S. ( 



subject to the available stowage constraint. It should be noted that 

each item could have different demand distributions and the model 

would still be valid. The only requirement for utilizing the model 

is that the density functions be continuous. * 

Table (1) contains the parameter values pertinent to the 

example. In accordance with the procedures of the previous section, 

an available tender capacity of 10,152 is assumed. It should be 

apparent that due to the iterative nature of the solution, 

will be approximately equal to 10,152 for an optimum set Z. 

(The most obvious reason is the need for rounding off the values 

obtained to intergers. In order to insure that the space constraint 

is met, the z^’s obtained should be rounded off to the lower interger 

Ni * 

with the result that ~ ). 

u: \ 

The values assigned to the items listed in Table (1) reflect 
varying combinations of essentiality, cube, expected demand and 



^Another model could be developed which would utilize discrete 
functions . 
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The essentiality valves in Table 1 were arbitrarily chosen to 
show that items 1-8 are approximately four times as essential as 
items 9-16. In a real life problem involving a large number nf 
items, a scheme using the military essentiality codes could be adapted 
for use in this model. 
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In order to determine the value of M once an optimum set Z 
is obtained, equation (3) must be utilized. For the case when 
demand is characterized by the normal distribution. 



2- oo 

(3) s t cz L l = -V- f.bOcLx; = 



12 




















' 





















O') 









-M^V 



Tzw 



-^zvr 



Appendix II contains a detailed derivation of equation ( 3 1 - ) . 

Since equation (5) requires z. 0 this implies that for the 

case of the normal distribution equation (5) and ( 3 1 ) are only valid 
. . 2 

when — > 3.87. (Fig. 1) 







Fig. 1: Graph of the normal distribution function 

Since => Zi ir^ >-3.8T. 

* For ‘Z.^> O , the ratio ^-> 3 87 must hold. Examination of 
Table 1 indicates that this condition has been satisfied. 

Due to the requirement above, the third term of equation (3 1 ) 



^ The value for-^- is dependent on the accuracy of the tables 
used. The value of 3.87 was obtained from the Chemical Rubber 
Company Standard Mathematical Tables. 
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can be ignored in this example since its contribution is insignificant 
for the values assumed. 



Table 2 is a presentation of the results obtained when the 
parameters shown in Table 1 are used in the model. The optimal 
set Z b =||^>,T-*. || which most nearly satisfies the arbitrary 

constraint of 10,152 cu. ft. tender capacity was found by using A~.04\6o 
Note that items 9-12 generally lag in percent of expected average 
demand filled &ytaa, for each value of A . This is intuitively 
appealing because these are items with low military essentiality and 
high cube. The items with high essentiality and low cube (5-8) are 
stocked at consistently high levels. 

Figures 2a thru 2d are plots of percent of expected average 
number of demands filled, versus a characteristic number P, defined 



as 



Mi 

C <Tc ' 



Since each item is characterized by the parameters m, <3“, M and c, 
each item is also characterized by P. Note that for different values 
of C 1 the curve shifts somewhat but retains the same general shape. 
In each case as the characteristic number increases, the percent of 



expected average number of demands filled increases. 

The effect of variability of demand for this example is investi- 
gated in Figure 3 by observing the behavior of items 9-12. Since 
items 9 and 10 (also 11 and 12) differ only in their values of 
we would expect any differences in their optimal stock levels ( Z-C ) 
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for various total cube requirements to be due to variability of demand 

(cr^ . 

Figure 3 supports the argument that when a large total cube 
is available, items with high values of <3" will be stocked in larger 
amounts than items with small <5~ . However, when the total cube 
available becomes more restricted, the stock of items with high 
variability of demand is reduced rapidly. Plots for other items which 
differ only in their value of G" show a similiar trend. 

Reference to Figure 4 indicates that as total available cube in- 
creases, the total essentiality of the load list increases also, but at 
a decreasing rate. Planners who wish to increase the spare parts 
capacity of a tender would have to weigh the cost of increasing 
capacity against the value of the added cargo. Use of a computerized 
version of the model with a large number of items would provide a 
new optimal load list for the proposed increased capacity. Reference 
to a curve similiar to that shown in Figure 4 would be an aid to the 
planner who must determine whether the essentiality of the added 
inventory is worth the cost of increasing tender spare parts capacity. 
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6. POLARIS Tender Re -supply (General) 

The maximum deterrent value of the FBM weapon system can 
be realized only if the maximum number of SSB(N)'s are maintained 
in the designated areas. The logistic support concept is based on 
satisfying this requirement. Lack of essential spare parts can be 
interpreted as a deficiency in the number of ready missiles scheduled 
to be on station. 

A load list is developed to allow the tender to provide effective 
support for a given period of time in the event re- supply is cut off 
for any reason. As a result, replenishment of the tender is in- 
fluenced by the following factors: 

(1) maintain stock levels at or near the optimum load list 
point , and 

(2) provide delivery of nonstandard demands (items not 
included on the load list for one reason or another) 
from the stock point to the tender. 

The factors considered in selecting the mode of transportation 
(air or surface) for an item include the transit time available , 
priority, physical characteristics of the material, capability and 
availability of the various carriers, security and cost. Of these 
factors those necessary to effect delivery to meet program require- 
ments should govern and considerations of cost and equitable traffic 
distribution should be regarded as secondary. 
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In order to illustrate the considerations mentioned above con- 



sider the case where a tender is assigned to support X SSB(N)'s. 
The load list for the tender is designed to provide Y days support 
during which time each SSB(N) is alongside for Z days. (Fig 5). 

Y can be interpreted as one cycle during which an SSB(N) has an 
upkeep period of Z days and a patrol of W days (Y = Z+W). If the 
assumption is made that the alongside availability period for each 
SSB(N) is the same, then 3 K is the number of SSBfN^s at 

the tender at any time- Under this assumption an SSB(N) departs 
and a new one arrives every days. When demand by each 



SSB(N) is assumed to be uniform, stock levels aboard the tender 



where there is no re- supply, the tender is completely depleted of 
spare parts at the end of Y days. 




Y 






Z. — H-* 



-w 



H 





Tender Support Cycle 
Fig. 5 



25 



7. 



Routine Re- supply 



Consider the first factor influencing re- supply mentioned 
previously. The average percentage of the optimum amount ('Z.'v. ) 
which the tender maintains will be defined as the effective inventory 
level. This level is a function of demand, re-order policy, and the 
lead time associated with the transportation of the item. Consider 
the case where demand for all items is uniform over a ninety-day 
period. Fig (6) is a plot of effective inventory levels at the tender 
versus lead time required for delivery as a function of re-order policy. 

The graph can be utilized in several different ways. If a re- 
order policy is chosen the graph indicates the inventory level that 
can be maintained for various items as a function of lead time. If 
lead times are fixed and a certain inventory level is desired to be 
maintained then the re-order policy can be obtained from the graph. 

It should be noted that due to the capacity of the ship and the 
lead time required for delivery of an item the tender will never 
average an amount equal to the load list quantity for that item. 

Several examples will serve to illustrate the use of Fig (6). 

(a) Consider the case where an effective inventory level 
of 90 per cent is desired for an item. The graph 
indicates that this level can be maintained by re- 
ordering when the stock level drops to 90% provided 
the lead time is 10 days. Another possibility is to 
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re-order when stock levels drop to 95% if lead time 
is 16 days. The re-order quantity in each case is the 
demand to date plus expected demand over the lead 
time . 

(b) If the lead time for delivery of an item is fixed then 

the re-order policy determines the effectiveness levels 
maintained. With a lead time of 12 days a 60% policy 
results in an effective level of 76%, an 80% policy 
results in a level of 84.5%, a 90% policy in an 89% 
level and a 95% re-order policy in a level of 93%. 

If the lead time is 30 days, which is analogous to that of surface 
re- supply by AK ^then the maximum possible effectiveness level is 
83% for the example shown. 

It is obvious that high inventory levels can be maintained 
effectively only for the case when the lead times are short. This 
illustrates the requirement for air transportation of spare parts as 
opposed to surface re-supply. 

In addition to the effect on effective inventory levels maintained, 
air re- supply has the advantage of reducing the confusion factor at 
the tender when supplies are received. Rather than being saturated 
with parts delivered by an AK all at one time the tender is able to 
maintain better control over spare parts storage. This in turn tends 
to reduce "loss" of items due to misplacement. 
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Lead Time 



8. 



Determination of resupply method for NIS/Nonstandard items. 
Consider a spare part that is required by an SSB(N) and that is 



not in stock (NIS) on the tender. If it is assumed that a spare part 
NIS on the tender is always available at the stock point then the 
mode of re- supply transportation (air or surface) can be determined 
by using a criterion developed as follows: 



Let , - 


Julian date of departure of AK from stock point. 


= 


Days AK spends in transit to tender. 


li 


Days remaining until j— SSB(N) starts patrol. 


U' * 
II 


Julian day that j— SSB(N) generates demand 
for an NIS spare part. 


c = 


Days required for loading, unloading and 
installation of spare part. 


LTT = 


Days required by spare part to travel by air 
from stock point to tender (MATS, commercial, 
etc). 



If — ^Av^-V-~^\vc -V- C* , re- supply the spare part by AK; 

If Rj -*-Sj < ^ak-V-Tak+C. > re-supply the spare part by air, 

subject to the constraint ' • 

If -v-S: < \_TT~ , then it is impossible to re-supply the spare 

J sj 

part prior to the SSB(N)'s scheduled departure. 
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9. 



Conclusions 



1. The number of stockouts will be reduced if stocks are 
maintained in accordance with the model. The model maximizes 
the average number of demands filled over some length of time. 
This is equivalent to minimizing the number of stockouts over the 
same length of time. 

Z. As the spare parts stowage capacity of the tender is in- 
creased, the total worth (military essentiality) of the inventory is 
increased, but at a decreasing rate. For example, if a small 
tender and a large tender each have their cubic capacity increased 
by the same constant amount, the small tender has a larger per- 
centage gain in the military essentiality of cargo than the large 
tender . 

3. Items with high essentiality to cube ratios are generally 
carried at high stock levels. For those items with nearly the same 
pattern of demand the model maintains high inventories of 

high essentiality, low cube items, and low inventories at low 
essentiality, high cube items. 

* 4. When items must be deleted from inventory because of 

space limitations, those with the lowest essentiality to cube ratio 
are the first to be eliminated. 

5. The spare parts stowage capacity of the tender is the 
limiting factor in determining the length of the period over which 
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demands filled should be maximized. The load list should be designed 
to maximize demands filled over a given length of time. 

6. Re -supply by air is required to maintain high inventory 
levels of essential spare parts aboard the tender. In order to main- 
tain essential spare parts inventories close to the optimum load list 
level at all times, it is necessary to have short re- supply lead times, 
or a large tender capacity. Since tender capacity is limited, high 
inventories can only be achieved by utilizing the short lead times 
associated with air re- supply. 

7. Re- supply ships cannot meet demands for non-standard 
items. The long lead time associated with ship re- supply precludes 
the possibility of meeting demands for non-standard items as they 
occur. For example, it is a rare occurrence when a submarine can 
execute an order for a non-standard item and have it delivered by 
ship re-supply before going on patrol again. 

8. The iterative nature of the model and the large number 

of items in a real life inventory make the use of a computer mandatory. 
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APPENDIX I 



Derivation of Optimal Tender Load List Model 




Essentiality of a spare part is an important factor in deciding 
on a load list which will maximize the average number of demands 
filled by the tender. On the basis of past usage data the model 
should indicate the quantity of each item to be stocked such that the 
total essentiality function M is a maximum subject to the tender 
stowage limitations. 

The function X^ ) will be introduced to represent a 

utility valve associated with the ability of the tender to meet demand 



It should be apparent that if demand is less than what is 
stocked the tender is able to provide 100% service. On the other 
hand if demand exceeds the stock level then the tender does as best 
as it can by providing an amount ~Z-C. • 



for an item. 



c 



ML Xc. X^^-2.; 

X(. >, "2. i 






L — M 
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In accordance with the measure of effectiveness mentioned 



previously the objective is to maximize the total utility value as 
defined above 

Max. 7. ...... Z. W V ElLj- 

subject to the constraint 

C ^"Z v ~ C- ■> ^ > O • 

^.Cx^-z^Y^ J JCx* > y i ..x N > )ci)( t ...c\x K i 

where \> ^31 ^ is the joint probability function of 

By assumption (7) this joint probablity density function can be 
written as: 

$ Cx v ,x x ^ f^Cx^ 

C<3 

L B v definition S 9 fwUC>A-xL =13 

As a result can be expressed as follows: 

ECz^uc.xcQ = y o ■ 

Let M; S.fc't = 

M; ^ ^ w Cx;,2^f v Cx^<Sx; 

Maximize ^ v. ^ v.C2. 

Subject to: - C_ 

L 

The problem presented is one of maximizing a function of v\. 

variables where the variables are related by a single constraint 
equation. The problem could possibly be solved by eliminating some 
of the variables using the constraint equation. Eventually the problem 
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could probably be reduced to an ordinary maximization problem. 
This procedure is not always feasible and therefore the technique 
of Lagrangian multipliers will be utilized. 

To find the points of where the 

variables "Z^Z*.. . -"Z.^ are related by the equation — C= O 

the function L is defined as follows: 

V - M 



where X is a Lagrangian multiplier 

O - XCc^-2,1- CH) 

= — Xcc , 

C v_= \, 2., M 

As a result^the following set of(n + Inequations in(n + 1^ 

variables ‘Z. v> Z i ...'ZNj and X is obtained. 



_ 







X (c.»-z. v - C_') - O 

o 






*4 



C. JZ; 






X ■$— ( cX - o 

- C - o 



By rearranging the terms the above equations can be put in the 



form 



( 6 ) 



jeL_ \Lmc. adz;. VI - Xcc 
dz; 
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Recalling equation (3) 

( 6 1 ) 

Lebnitz's rule is utilized to evaluate the left hand side of this 
equation since the units of integration depend on the parameter 'Zi . 

Rule : V WCXjVjdfc = > tl b'(V) - V\_Ca-CO , 

-4- r cL* 

Jg-CM 

^_C5 o X.f v U^ckx w 4-2c^ ^ 

= C^w^u^-cm-o] + Co- 

= L. £l<^0<kx^ 

Since r frltxC) &-xi ~ IL - X 0 .fcCxO &xl -t- j* 2 , ^^CxOcCxc. 

C*^5 "2. * 

Then J £(jO&XC. - 1 — J“ iLOO&XC 

Z «L o 

By substituting and rearranging terms equation (6 1 ) reduces 
to the following: 

1 — S Q 5^bO<$b<v - 

or J o = 1-XC^.) (5) 

Equation (5), subject to the constraint equation, provides the 
method where-by the set Z can be determined. The two equations 
determine a unique set of *Z^'s provided that for all L ^ ^CX^>Q 
for all values of . 

The distribution function , , of a numerical valued 

random phenomenon is defined as having as its value, at any real 
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number “Z_ 



, the probability that an observed value of the random 



phenomenon will be less than or equal to the number 2L . 

- PLZ- £ 2L~^ . Since the left side of (5) is a distribution function 
the following condition applies: 

o- ■C'fct-'Oi'ioc 4 1 . 




In order to solve for 21 L in equation (5) we must assume a 
value for A . It would materially decrease the iterative require- 
ments if the possible values of \ could be confined to a portion 



of the real line. 






Since 


o S*' JcUOAx;. & 


1 


Then 


° < i.- X( ^ 








> X 










-X(^< o =« 


> X ^ O 


• 

% • 


> A > o 





These limits apply for all v_ and therefore \ must fall in 
the range . ^ O in order for all items to be 

included in the load list. 
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APPENDIX II 



Derivation of S^C ’Zl') for the case where demand follows the normal 



distribution: 



where 

oo 



but 



S 0 5aO<£*-Xv - ~ S 0 £(xMki fcXxO 

.\ £ 1 - J*" f^GO^XC. 

4 v. O 

fcUOAocC. - 1- XC^ 

■'* ^ f.Cxa&*l «0~~ Cl-A^J}} 

.% ^Cz;^ = X^xcfAxo^*; -v -z^Lv-cv-a-^q 

In order to solve for an explicit expression must be 

found for V*iSc Cx^^lxc. . The procedure is as follows: 

■z t x — -1 C 

(II- 1) S a Xi&LXl)&-XL ~ S Q ^ " C ^ XC - 



± -1l ^ 

By taking the derivative of -==. S_ 



X. - Vvx„ 



■) 



\7^r 



the 



following expression is obtained: 



a J_ ~ 

-^TV^" 



^ ^ V cr v n^Z 



L / X v - vv\ v \ x 

* v. -5=r ) 



- _(,Xv.-wO ^ 

~j-tt <3V*- 

- _ x: e ~ ^ ^ -L- O ^ <r ^' 

-?^R X -^rcr; 

Integrating both sides, of the equation between the limits of O 

4 

and ~Z> ^ : 



-,^-W 

c a _i_ ( 

J o 3x;^r 



>■ / XC-wva ^ x; 
- x v <3-^ 1 r 



x/ x:- K .V 
- A o-v. >. 






r : -tC^V r 2 " . 

J 0 ~^TT^ 1 
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5 



Rearranging terms: 
_x( x v -m-y 






-\xrr<5V 



d*. 



•z-c 

= 5 , 



-tc^v 

«• ^ 



v>U 



-U 






c3o< 



Note that the expression on the left is the same as equation 
II- 1 except for a factor OT^ . In order to obtain the same ex- 
pression multiply both sides by OT^ which results in 



s 












•A^ivCT; 



-TV 

6 



y ( x^ :r z ' 

rt» r- - . - ^ \~5T~ 1 r ^ -xV cr. » 

^-=X ^r 2 - *■* ■ w- 



<=1 



* ( 

_ _xV— » 

but _ \ d Sfe e- = j=- 

J 0 dx* -vxrr 

-i- f-zc-^OX* 

- 5=r ' 



<S~c 

-\"wr 



-tc^r 



— *_ v cr:- * 



G- 



G~C — z_ v <rc 
' VXTT 



Observe that 



v ( X^-W\ A -2C 

-^v cr, r \ 

S 0 jW'T-Z- 



A. ( X v-^ V~ 

_ V <3- ' 



but 



X" fctxJ<W = t \ 

.*. Vv\ 4 C~ '-4- e,” ^ ^ ^Xv = K vC\-X(~;^ 

d Xlxrcr, 



Combining the above derived relationships an explicit equation for 
results : 



(3 .) s,t*.L^ = ^cD - ^ x 4r i ''i 

.a; ^ 

for "Z.^ >0 



38 



APPENDIX III 



A computer program based on the flow diagram illustrated in 
Figure 7 was written in the Neliac language and compiled on the 
United States Naval Postgraduate School CDC 1604 digital computer. 
The Neliac statements are found in this appendix following Figure 7. ^ 
Outputs for several choices of parameter values confirmed the re- 
sults shown inS ection 5. 

As indicated by the computer flow diagram in Figure 7 the 
first step in the operation of the program is to carry out the routine 
input procedures. 

The next step is the computation of the value of lambda equal 

/ \ 

to the smallest fraction, v ) . Denote this minimum ratio as 

min ^ -q~ / 



This value of lambda is used to compute all values 



\ C 

of AREA = ( ,i_ — X -pjr ) for i = 1, . . .N. Each ~Z.[_ is computed 



from: 



'vjvi. - 



-^r ( =f£') X ' 



which provides an approximation of the normal distribution function. 
In the expression above: C - X. ^>8 

where " \ -V O- \?357L\ ^ 

-V- C.02J1H3^VS (^-O.C003 c \H%(^] 

-v O. * 



^Details of the Neliac language may be obtained from Prof. 
R. M. Thatcher, Operations Research Department, U.S. Naval 
Postgraduate School, Monterey, California. 



39 




















' 





























































The product C^.'Z.C. is computed, stored and summed. The result, 
<^L~Z-C is stored in the cell named SUMCZ. 

This sum, , is then compared with the contents of 

Big C , the specified total cube available for the storage of spare 
parts. If SUMCZ > Big C, then the cube of the load list, for the 
given lambda, is greater than the capacity allowed for storage. When 
this occurs the item with the lowest essentiality to cube ratio is dis- 
carded, a new X = minl^c^ computed from the remaining items 
and is used to determine a new SUMCZ. This new value of SUMCZ 
is compared with Big C. If SUMCZ remains larger than Big C the 
item with the lowest essentiality to cube ratio is discarded. This 
procedure is repeated until SUMCZ becomes equal to or less than 
the cube available, Big C. When a new load list is found that meets 
the cube constraint then the essentiality of the load list is computed 
and a print -out is obtained for Z; ^ A ^ an d 

A table of the items discarded is also provided. 

If the cube of the load list computed for the initial value of 
lambda ( vvou^ ) is less than the capacity allowed for 

storage, then a print -out is obtained for an d 

In order to find the load list that most nearly meets the cube constraint 
it is necessary to evaluate for several assumed values 

of X 
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The number of values of A to be assumed is specified in the 
program by the parameter K. If, for example, K~20 then the 
variables and will be evaluated for 

W V 

\ w ' 

twenty equally spaced values of A between min ^j -r- and zero. 

In this case twenty print-outs will be provided. The optimum load list 
is taken as the load list that most nearly satisfies the cube constraint. 
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Fig 1 , 



Computer Flow Diagram 
42 
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Neliac Statements for Computer Program 



= vt, 

= 






( 00.0 x o 


(00.0 XO- 


S’O.O X o , 


SO.O^O, 


Voo O X o , 


(oo.o x o, 


50.0 X O, 


so. oxo, 


\co.o )to, 


(00.0 X o , 


50-0 xo , 


50-0*0, 


\oo .0X0, 


VOO. o x o , 


so.o xo , 


C05Q.O XO , 


\7_.0 x O , 


3-0 X O- 


WO x o , 
V3.0 x o , 


3.0 x o , 

3-0 x O, 


V2..0 X O , 


3.0 X o , 


VL.O X O 5 


3.0 X o , 


\7L. O * O , 


30 X O, 


W. O X o , 


30 XO, 


W.O <o , 


0303. 0 *-0 , 



VJ (U^ - A-O * O .. AO * O . 



A 



COb^i - 



Ao * o j 

Cy.O * O ^ 
A-O KO , 
Z.O * O , 

z.o *o , 
z.o *o, 
•Z..O *o , 

VJL.O * O 

\5L.o * o 



AO * O, 

AO X O .. 

Ao * o’ 

"5-0 * o, 

" 2_-0 * O , 

7 - O * O . 
COoTL.O VC O 



3-o 

"5.0 

V2_.0 

VT..O 



X o 

X O 

* o 

* o 



3.0 * o 

3.0 * o 



A 



) 

•} 

3 

3 

> 

•) 

3 



VZ..O XO , 
WO * O, 

3.0 * O . 

3.0 **c o , 
vx. o * o , 

V2.. O *0 , 

3. 0 * O , 
Oco'i.Q VC O 



xrej-\ wtvt'N - Vj z., s, a,s, <o,n, s, s, 

\0,V\,«., V3, VA , VS-, VC^ 

axk\ z o x o 

\<K - ■20. 0 X O 

XX- OQO, 

LAHSQC\ = CCOO.OCcC *. O , 

Sonc.^. - oooooo . ccoo x o , 

SOK\\kJS^ 2_ “ OOOOOOcOOCO ^ O J 

AR.eA = O X O J 

"2_ (v<=^ r oooooo. ccoo X O , 

S^E-Cvu,^ - Oooooo - Oooo X O, 



/\ Indicates a blank card 
3 Indicates a card of all nines 
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A 



A 

/\ 



A 



A 



A 

A 



A 



VJOC- ~ 0*0 , ^REAn^ = Ooo.OoaCi^O, 

\GOGaTT ~ — . S' * o 1 

K~ 0*0, Y ~ o * O , T" = ooo. oooo * O . 

svGv 4) zj>r —o*o 

E.KP1 ■= O * O , 

EXPZ -0*0, 

SQRT^px = 0*0 

TvJc|)pT = b^ssi^slonz.xo, "XX^OQ, 

FAFv_AQ J K\G\ > Mvvsi , 

Fv_AS, ^ TF\J\<Ei j XTt_k\ YnsCARO AO| ~ Ooo , 

~rL =0*0 . tx. -oxo, 

*4oux_ - i~m ~? oooTrrn ooo & , 

; START” i , {;<>-, \ , 

, 1 « Pfctfi<5RAR\ t _ OArA\CtTb\ Pa^E-X » , , , i 
PRIMT- PARWAj ’ 

O — ^ VTAC* — *- TFYaG^ — »~ WFv./^ — K\ 

Sqft,TF ( T\.>i<JiF»TL •, ^CSi'CtTZ.TTA , 

F\KiO Hiv4 '•A’-qbC. ; 

o — *, VA , 

^£P€AT ; 

TF m~X\ MR O'A^J - FAOV_\_ FA -V 1. — »• v\ } v\ = Vvitvi : 
_ Anush . tropokt-. •* 

\MOa 1 / c. Cnj *- LAMSCq , 

O 2- VAXM ''Uv5\ , 

F<jbR. K - Otq vAkn — i_ 

"XP VTVOA VM«. OG - •. Ov'SCA.Q.O Vsitpc. . 

\M VAk.A / — a- \m^c < : 

l -AM\ 50 A V< — *■ VAVMi F\AG\ * ) 

OvSCAFO VAC^C. V > 



IF MFLA.4 = O ■ UASOA 



-*■ MIM 






> •) 



tOAPOTvS. ArssA, : 

o -■> ■S'iU'AC. 2 -. , 

F#>R 3T - Olx\ - 1_ 

(. ‘XF 3ZF3A MHCXl - Mi Ou_ 






i-Av\QO A * C.CCX3 /xur.TG^ — ^ AftOA CX^i 
x — rA - <~s • o— 2, . Q~ -i— 1 



O-S 



- Ars~A 'CxA 
c-^tAV. r 



1*. O • 

TF A CXT = O *. 

T LXG ^ O-S 

AROAt-xS <- o.ocsoojl : O — *-A 
wormeX (rrsam5 \TA 
O — T — *“ TX C-^nTT. 

„ „ _ A P.EACXl - O- F — ARtAti 3 , 

APvM-<^v\j 5 - <_ o. ooool *. O — *» 'X, cowrr. ; 

HORHEX, ( . AROAVA5 ; 'T A * 



OfcMTT. 



■AuOxh' 



Ojfc»F4T I 

X "V \ '* n -* , ** i 5^ ^ ^ 

J i « TX \ - V ^ XX A VJ UMStlA \ * \ 

> l_AvAfo\j/\ 4 \j A.CEA \ - V > AftRA XXI < 
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\J AREALS \ 

= \ > 5 VJ <~T \ = \ >1“ > 1 

KKCXl -v- SVGVAA * T~ ^ *S- CX3 

2.CX1 * c.Cx.“X -v ^'o^e=a_ — *» 2 vjvac_"Zl , 
<£va\t *. ^ , 

, ^■prvx'TT « Sut^ct \ — \ > sonnc^sl > ^ 

32F TFv_A<^ — O •- 1^.5T, c.^MPuTS - 



A 

A 



A 



A 

A 

A 



A 



VEST-; ^xotAC_r2; <- B\GkC_ *• C^nVuT^S^. j 

1 — - FVAQ , 

K>G* -V- 1. »> 5 

(A3TSA F\_AGv =»r \ — » VTTXM X>\SCAOJ^ X ''■O. , 

VO -V i_ — *. VO , 

NiUUV_ *- X\BA VOVR X.V_“i j 

f\^0 VAXVsi VOcfc>C_. 



a— ^ nr 

Q.^>VAPOT-<El 



's*fc *. , <> , > ? > , 

Ua\x VJVJ "2: OvJ S » , } 



O 



Sonwl , 



) 



<Jjvo : f<£>r. X - o(.i^ NiNi — i_ {_ xf \TeitA voRtixl - vauvao 5<\p, ; 

“SvC^mAvXxOl / SQRT 2. RE — S\Qc^2_PX , 



VAK.C.X 1 / 5 \<aRA CX3 — ■— X X x ®- X - 

RES^S X X — ► X , BXPF (X EXP1), 

(I.CT3- vaxcx! /s\gvaaCx3 — »■ Y* Y— »-Y, 
vavxgrt S" * X expr (A -, expz'1 , 

m\cx3 * arcaCx! vz.cxn * U-o -A rea tso') 
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